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1. Introduction 

The lattice study of pure-glue gauge theories gives interesting information on the behavior of 
non-Abelian gauge theories, at a fraction of the typical computational costs of lattice QCD, and 
offers the possibility to compare the results of Monte Carlo numerical integrations with analytical 
calculations. Here we present a summary of our recent work [1] in which, following this line of 
research, we studied the equilibrium thermodynamics of SU(2) pure-glue theory in its confining 
phase, i.e. for temperatures T below the critical deconfinement temperature Tc. We show that 
the equation of state in the T < Tc region can be modeled as a gas of non-interacting relativistic 
glueballs, provided that the contribution of heavier glueball states is described in terms of a bosonic 
string model. This work can be considered as a generalization of the study presented in ref. [2] 
for SU(3) Yang-Mills theory, later extended to SU(A) theories in 2-1-1 spacetime dimensions [3]. 
Related ideas have also been discussed in refs. [4, 5]. The motivation to focus on the theory with 
N = 2 color charges is that it provides a crucial cross-check for the string model, as it admits only 
states with charge conjugation C = -|-1 and it is characterized by a second-order deconfinement 
transition. Furthermore we show that similarly good agreement (with no free parameters) is also 
found for the SU(3) theory, using lattice data computed in ref. [6]. 


2. Thermodynamics on the lattice 


A quantity of major phenomenological interest in finite-temperature field theory is the pressure 
p, which in the thermodynamic limit F —)• oo equals the opposite of the free-energy density/: 


p = — lim f = lim — InZ. 

V—>oo V— V 


( 2 . 1 ) 


The pressure p is related to A, the trace of the energy-momentum tensor (also called trace anomaly) 


j4 gp \ p4 


( 2 . 2 ) 


Two other quantities such as energy density and entropy density can be readily evaluated: 

* e InZ A-|-4n 

= A-|-3p, ‘S'= — H—-—= ———. (2.3) 


dlnZ 

e = — 


V dT 


The SU(2) Yang-Mills gauge theory is regularized on a four-dimensional hypercubic lattice A of 
spacing a and hypervolume a^{Nl x Nt), with periodic boundary conditions on all directions. The 
temperature of the theory, according to thermal field theory, is the inverse of the shortest (temporal) 
size, i.e. T = \/{aNt). In this work, variations of the temperature are performed by changing the 
lattice spacing a (which is a function of the coupling) while keeping Nt fixed. The lattice version 
of the action is set to be the standard Wilson action 

Sw = -4L L Trf/^vW (2.4) 

^ xeA0<li<v<3 

where g is the bare lattice coupling (which is related to the Wilson parameter j3 = 4/g^) and U^y{x) 
denotes the plaquette from the site x in the (ft, v) plane. 
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Figure 1: Values of string tension in units of a obtained from our lattice simulations are showed along with 
those reported in ref. [10]. The solid black curve shows the interpolation to the functional form in eq. (3.2) 
with the associated uncertainties. 


The dynamics of the lattice system is defined by the partition function 

x6A/i=0 

SO that the expectation value of a generic, gauge-invariant quantity A is given by 

^ xeA/i=o 


(2.5) 


(2.6) 


Any expectation value is estimated numerically via Monte Carlo numerical integration averaging 
on a large set of configurations generated by a mix of “heat-bath” and “overrelaxation” algorithms. 

Thermodynamic quantities can be obtained from plaquette expectation values by the “integral 
method” [7]: the pressure (with respect to the value it takes at T = 0) is given by 

p = ^ [\i5>{{U^)t-{U,)o) (2.7) 

r Jo 

where (t/p)r denotes the average plaquette at a generic temperature T. The integrand in eq. (2.7) is 
closely related to the trace anomaly A, since: 


A = 


6 dp 

a'* d\na 


((t/p)o- 


mT), 


(2.8) 


where /difna) can be readily evaluated from the scale setting. 


3. Scale setting 

In order to determine the temperature reliably at a certain value of j3, a precise scale setting of 
the theory is mandatory. We computed non-perturbatively the zero-temperature Polyakov loop cor¬ 
relation function, denoted as G{r), for different values of r and j3, using the multilevel algorithm [8] 
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on 32^ lattices. The interquark potential T(r) was extracted from T(r) = —[lnG(r)]/(aA^,) and fit¬ 
ted to the functional form 

Ka 

aV = aor + uVq — (3.1) 

12r 

using the tree-level improved definition of the distance r [9] to obtain the string tension a in units 
of the squared inverse lattice spacing. As a final sfep, we performed a polynomial inferpolafion for 
fhe logarifhm of fhe siring tension for differenl values of fhe Wilson parameter j3 

ttpar^ 1 

log(aa^)= Yj (3.2) 

AO 

wifh = 4 and jSq = 2.35. The resull, which models fhe relafion belween a and j3, is shown in 
fig. 1, and allows an accurate determination of fhe femperafure for a large range of j3. 


4. Numerical results and comparison with a bosonic string model 


The main parf of our numerical sfudy of fhe SU(2) Yang-Mills fheory is focused on fhe equa¬ 
tion of sfafe in fhe confining phase of fhe fheory. The resulfs for A/T^ were obfained via Monte 
Carlo simulafions on lattices wifh differenl temporal exlenls, keeping fhe aspecf ratio Ns/Nt large 
enough lo avoid finile-volume effecls. The resulfs are showed in figure 2 and are plotted againsf 
T/Tc using Tcjs/o = 0.7091(36) from ref. [10]. The only physical degrees of freedom of fhe 
fheory in fhe T<T, region are massive glueballs: if is reasonable lo assume lhal such slales are 
weakly inleracling wifh each olher' and Ihus fhe system can be modelled wifh good approximation 
as a free, relalivislic Bose gas. The Irace anomaly of fhe latter is given by 


m^T “ Ki (nm/T) 

A = —L- 


271 ^ 


n—[ 


and using asymplolic expressions for fhe modified Bessel funclion we have 

3T 


, , exp(—nm/T) 

A ~ m I — 1 Y - 

n—\ 


271 ) 


,3/2 


1 + 




(4.1) 


(4.2) 


In figure 2 fhe conlribulions due lo fhe lighlesl glueball only (fhe 0++ sfafe) and fo all fhe 
slales wifh mass m < 2mo++ (faken from fhe speclrum calculaled in ref. [13]) are shown. The mosl 
sinking fealure of fhe figure is fhe large mismafch belween fhe glueball gas predicfion and lattice 
dafa for T close fo Tc. To address Ibis mismafch, we assume lhal fhe densify of heavier glueball 
slales is described by a Hagedorn speclrum: indeed, only an exponenlially increasing speclrum can 
accounf for fhe exponenfial dependence in eq. (4.2). In particular, a Hagedorn-like speclrum arises 
if we model glueball slales as Ihin closed color flux lubes lhal can be described in ferms of closed 
bosonic slrings [14, 15]. Specifically, such closed bosonic siring model leads fo a speclral densify 
(see ref. [3, appendix] for a derivalion) 

_exp(m/T») (43) 

*The expectation that glueballs are weakly interacting is borne out of theoretical arguments in the large-A limit, but 
lattice results indicate that these expectations are surprisingly accurate even for the theories with N = 3 or N = 2 color 
charges [11, 12]. 
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Figure 2: Comparison between our lattice results for the trace anomaly in SU(2) Yang-Mills theory from 
simulations with different Nt and the behavior expected for a gas of free, massive glueballs. The dotted 
line cotTesponds to the contribution of the lightest state only, with quantum numbers = 0++, while the 
dashed line includes all the low lying glueballs with masses lower than 2mQ++. The solid line includes also 
the contribution from high-lying states, described by a bosonic string model. 


where the only free parameter is the Hagedorn temperature Th [16]. If the effective action governing 
the string model is identified with the Nambu-Gotd action [17, 18], then is fixed and its value is 


7h = Tng = 



0.691v^- 


(4.4) 


For SU(2) Yang-Mills theory, however, the deconfinement transition is second order and the Hage- 
dom temperature coincides with Tc, so that no determination of T„ is required. Whether the transi¬ 
tion is continuous or not, the contribution of the complete glueball spectrum for a thermodynamic 
quantity such as the trace anomaly A can be written as 

noo 

M'^)= (27-1-l)A(m,-,r)+nc / dm'p(m')A(m',r), (4.5) 


where the first term includes the contribution of low-lying states (whose masses are taken from 
independent lattice calculations) up to a threshold chosen as rrii^ = 2mo++, and the second term ap¬ 
proximates the contribution of heavier states via the bosonic string spectral density The (pseudo)real 
nature of the representations of the SU(2) Lie group allows only glueball states with quantum num¬ 
ber C = 4-1: thus the multiplicity factor nc is set to be 1. The final result can be seen in figure 2: 
lattice data and the bosonic string model prediction are in remarkable agreement. 

We tested this model also for the SU(3) theory, using the results for the equation of state from 
ref. [6] and the glueball spectrum from ref. [19]. The differences from the N = 2 case are: 


• SU(3) Yang-Mills theory allows for both C = ±1 states, thus the factor nc in eq. (4.5) is set 
to 2 in order to account for this charge-conjugation multiplicity; 
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Figure 3; Comparison between the prediction of a massive-glueball gas, including the contribution from 
states modelled by a closed Nambu-Goto string model, like in eq. (4.5), and continuum-extrapolated data 
obtained in ref. [6] for SU(3) Yang-Mills theory, as a function of T/Tn. Our data and predictions for the 
SU(2) theory are also showed for comparison. 


• the deconfinement transition is of first order and the value of the Hagedorn temperature is 
fixed by the Nambu-Goto prediction, see eq. (4.4). 

The resulting curve for the trace is showed in figure 3 as a function of T/Tn, along with the results 
for SU(2). The parameter Tcjs/o = 0.629(3) from ref. [20] was used to plot the data; this value 
is consistent with the recent result roT^ = 0.7457(45) taken from ref. [21] and combined with 
ros/o = 1.192(10) from ref. [22]. Excellent agreement with the effective string prediction is also 
found for SU(3) Yang-Mills theory, and the contribution of C = — 1 states in the spectrum is crucial. 
We remark that the non-interacting glueball gas predictions for both N = 2 and N = 3 colors do 
not depend on any free parameters: furthermore, these findings are in agreement with previous 
results in D = 2 + 1 spacetime dimensions for SU(N = 2,..., 6) theories [3] and can be considered 
a generalization of those. 
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